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f->, Abstract: When solving renormalisation group equations in a quantum field theory, one 

^ often specifies the boundary conditions at multiple renormalisation scales, such as the weak 

I— I and grand-unified scales in a theory beyond the standard model. A point in the parameter 

,__! space of such a model is usually specified by the values of couplings at these boundaries of 

^ the renormalisation group flow, but there is no theorem guaranteeing that such a point has 

\^ a unique solution to the associated differential equations, and so there may exist multiple, 

■^ phenomenologically distinct solutions, all corresponding to the same point in parameter 

space. We show that this is indeed the case in the constrained minimal supersymmetric 

^s standard model (CMSSM), and we exhibit such solutions, which cannot be obtained using 

cn out-of-the-box computer programs in the public domain. Some of the multiple solutions we 

exhibit have CP-even lightest Higgs mass predictions between 124 and 126 GeV. Without an 

exhaustive 8-dimensional MSSM parameter scan per CMSSM parameter point to capture 

^\ all of the multiple solutions, CMSSM phenomenological analyses are incomplete. 
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1 Introduction 

It is a truth universally acknowledged, that renormalisation group (RG) flows are unique, 
once a boundary condition for each coupling involved in the flow has been specified. 

Like many universally acknowledged truths, this one is not necessarily acknowledged 
in the Universe inhabited by mathematicians. Indeed, the closest one can get in terms 
of a theorem (called Cauchy-Lipschitz by francophones, but due to Picard-Lindelof [1]) 
concerns the uniqueness (and existence) of the solution to the initial value problem of a 
sufficiently well-behaved system of differential equations in a neighbourhood of the starting 
point. There are many situations in physics where these conditions are not satisfied and 
so the issue of non-uniqueness (as well as non-existence) rears its ugly head. 

One physical situation where non-uniqueness is manifest arises in Sturm-Liouville prob- 
lems, namely linear, second-order, ordinary differential equations (ODEs) with homoge- 
neous boundary conditions specified on either side of an interval. These may be thought 
of as a pair of first-order ODEs and, as every undergraduate knows, there exists either 1 
(namely 0) or infinitely many solutions (namely an eigenfunction multiplied by an arbitrary 
constant) depending on whether the ODEs correspond to an eigenvalue or not. 

A different example, involving RG flows, is pertinent right now at the Large Hadron 
Collider (LHC). Disciples of high-scale supersymmetric models (rapidly becoming an en- 
dangered species) wish to know whether their models are ruled out or not. These models 
typically impose a large degree of uniflcation of the parameters of e.g. the Minimal Su- 
persymmetric Standard Model (MSSM) at a high energy scale, be they gauge couplings, 
soft supersymmetry (SUSY) breaking masses, or whatever. Such constraints play the role 
of boundary conditions at one end of the RG flow. The MSSM is then further constrained 
at the weak scale where various Standard Model measurements, such as the mass of the Z 



boson, also play the role of boundary conditions. The RG equations (RGEs) are non-linear 
ODEs, for which a solution satisfying both sets of boundary conditions is to be found. 

For concreteness, let us imagine an n-parameter flow (with n of O(IO^) in the MSSM), 
in which k measurements are performed at low energy, and there are / unification conditions 
at high energy. One may then attempt to answer the question of whether the theory is ruled 
out or not in the following way: choose values oin — k — l of the unified parameters at high 
energy, and call those values a 'point in the parameter space of the model'. Given that n 
boundary conditions have now been specified, and invoking the aforementioned universally- 
acknowledged truth that the flow is unique, one may find 'the' flow that satisfies the given 
boundary conditions, by means of a numerical iterative algorithm.^ The resulting flow then 
predicts the values of all particle masses and couplings of the theory, including those which 
have yet to be measured, but are subject to limits {e.g. limits on superpartner masses 
from LEP and LHC). If the limits are violated, the point in parameter space is ruled out. 
Finally, one can scan over points in parameter space ad nauseam. 

The problem with this approach is that no theorem guarantees that a solution found 
in this way is unique. (Indeed, the similarity with a Sturm-Liouville problem suggests 
the contrary.) Thus, there may be more than one trajectory satisfying the boundary 
conditions, each of which reproduces the unification conditions and the measured Standard 
Model parameters, but each of which may have completely different values of the as-yet 
unmeasured low-energy parameters. Some of these may be ruled out and some may not. 

This raises the spectre (horrifying or enchanting depending on one's spiritual taste) 
that points in the parameter space of, e.g., the CMSSM that have previously been ruled 
out, are in fact not ruled out at all, because there are multiple RG flows that correspond 
to the same point in CMSSM input parameter space, with one or more flows still allowed 
and yet not found by existing algorithms. 

In what follows, we shall find multiple RG flows corresponding to the same CMSSM 
input parameter point. ^ Though our numerical calculations hint that such points may be 
uncommon, it will become clear that we have no way of exhaustively categorising them. 
This, unfortunately, makes it rather difficult to say whether some regions of some high 
scale SUSY breaking scheme, e.g., the CMSSM, are ruled out or not. 

Before doing all this, we first try to convince the reader, in § 2, of the existence of 
multiple solutions for a rather more simple RG flow, namely the flow in the neighbourhood 
of the Berezinski-Kosterlitz-Thouless (BKT) phase transition. Then, in § 3, we exhibit and 
investigate the phenomenon for the CMSSM. It is not our purpose here to perform any 
detailed phenomenology. Rather, we content ourselves with pointing out the existence and 



^In a nutshell, such an algorithm works by choosing initial guesses for the a priori unknown parameters 
at, say, low energy in order to create an artificial initial value problem, which can be solved to find high 
energy values of all n parameters. Those which are known at high energy are discarded in favour of their 
known values, and one flows back again. The process is then repeated indefinitely, in the hope of converging 
on a solution. 

^Recently, Ref. [2] investigated different 'branches' of electroweak symmetry breaking in the CMSSM 
(and other models) given current constraints. However, these are not different multiple solution branches: 
for a given CMSSM point, only one solution was found. 



prevalence of multiple solutions in the CMSSM parameter space and investigating a few of 
their properties. We discuss the implications and context of our results in § 4. 

2 Toy example: The BKT phase transition 

To find a tractable example of the phenomenon of multiple solutions, we consider a system 
of RGEs with just two couplings, which is the minimal case in which the notion arises of 
imposing boundary conditions at different scales. 

A suitable example is the RG flow corresponding to the BKT phase transition [3, 4]. 
This phase transition was first studied in the context of the 2-d XY model in condensed 
matter physics and can be used to describe superfluid films and arrays of Josephson junc- 
tions. It also appears in particle physics in the deconfinement phase transition in compact 
C/(l) gauge theory in 2 space dimensions, where Polyakov showed long ago [5] that vortex 
configurations of the gauge field lead to confinement at zero temperature. The finite- 
temperature deconfinement phase transition was shown to be of BKT type by embedding 
in the Georgi-Glashow model [6, 7] and also by a variational analysis [8]. 

The BKT phase transition is rather special in that it exhibits a line of fixed points 
ending in a critical point; in the neighbourhood of this point, the RGEs may be written in 
the form 

dx 2 

1=... (2.1) 

These equations may be solved easily enough by a mathematical beginner, but physi- 
cists who know more than is good for them may find it amusing to convert the problem 
to a mechanical one, in the following way. The RGEs do not form a Hamiltonian system 
(since dxlp' 7^ —dyxy), but may be turned into one by the non-canonical transformation 
y = e^. Then z = x and x = e^^, such that we may think of a particle of unit mass and 
momentum x moving in the potential — ^e^^. The Lagrangian is t-independent and the 
conserved energy is, a la Nother, E = gl^;'^ — e"^^) = ^(^^ ~ v'^)- ^^ other words, the tra- 
jectories are hyperbolae in the (x, y) plane. Moreover, t-translation invariance implies that 
the conditions for the existence of solutions to the boundary value problem can depend 
only on the difference, ti — to, between the initial and final times. 

There is, in addition, a symmetry of the equations of motion that is not a symmetry 
of the Lagrangian, and therefore does not imply a conserved charge. This symmetry is the 
rescaling {x, z, t) — )■ (ax, z + a, a~^t). This symmetry implies that we can, without loss of 
generality, always rescale a finite time interval to be unity, ti — to = 1- 

The form of the potential, — ge^^, indicates that the momentum, x, must increase 
monotonically with t. There are thus 3 types of trajectory, as follows. Trajectories with 
x(to) < and positive energy, x^ — y^ > 0, proceed inevitably to 2; — )• — cxd, i.e. y = 0, 
in time — f °° / ^ ^z ~^ ^^- Trajectories with x(to) < and x^ — y^ < eventually 
turn around and proceed to z — )• +00 in time /°° , ^ ^^ < c«, as do trajectories with 
x^to) > 0. The trajectories are thus as shown in grey in Fig. 1. 
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Figure 1. Phase diagram and RG trajectories for the BKT system. The equations enjoy a scale 
invariance, so the figure has the same character at all scales. The phase diagram is indicated by the 
background flow-lines (in grey). We show (in black) two distinct sets of flows, with each beginning 
at a fixed value of x on a blue line, and ending after a fixed interval of time, At, on a red line. 
For the flows ending on the curved red line, At is finite, and multiple solutions are indicated by 
the intersection of this red line with the straight, dashed green line. For the flows ending on the 
horizontal red line, At is inflnite. 



It is apparent from the above discussion that the initial value problem is well defined 
for any finite starting point (x(to),y(io)) and finite time interval, but we wish to consider 
instead the mixed boundary value problem with, say, x(to) = —o, and y(ti) = b. It is 
immediately evident that the existence of a line of fixed points permits infinitely many 
solutions to the boundary value problem with 6 = on an infinite interval: all trajectories 
with x(to) = — a,y(io) < ^ will arrive at the fixed line y{ti — )• oo) = 0. This situation is 
depicted in Fig. 1 by the mapping of the smaller blue line to the horizontal red line. 

We are more interested in multiple solutions on a finite time interval. One can show 
(by straightforward, but tedious, consideration of the explicit solutions to Eq. (2.1)) that 
these never arise for boundary conditions of the specific form x(to) = —o. and y(ti) = b. To 
find multiple solutions on a finite time interval, we need to consider more general boundary 
conditions. Consider, for example, boundary conditions in which we fix x at to, and fix 
some (inhomogeneous) linear combination of x and y at time ti, say Qx(ti) + /3y(ti) = c. 
It is then easy to see that multiple solutions must arise for particular choices of a, /3 and 
c. Indeed, consider all flows that satisfy x(to) = —o-, for some other value of a, as per the 
longer blue line in Fig 1. If the flow were linear, the straight line x{tQ) = —a would be 
mapped into another straight line at t = ti, which would intersect once with the linear 
combination appearing in the boundary condition. (In special cases such as homogeneous 
boundary conditions, a = c = 0, there will be 1 or infinitely many intersections, as for a 
Sturm-Liouville problem.) But since the flow is non-linear, the line x(to) = — o is mapped 



to a curve (shown in red) , which will intersect multiple times with a suitably chosen straight 
line (shown in green) corresponding to the boundary condition at t = ti. Since, we have 
effectively converted the boundary value problem into an initial value problem for each 
point on the blue line, for which the RG flow is known to be unique, the number of solutions 
to the original boundary value problem is simply given by the number of intersections of 
the red curve with the green line. 

Notice how, in this way, we convert the problem of finding multiple solutions of a 
system of differential equations (the RGEs) into a problem of finding multiple solutions of 
a system of algebraic equations (namely the intersections of the line and curve in Fig.l). 
This makes it relatively straightforward to establish the existence of multiple solutions for 
a more general system (though we still need to be able to solve the non-linear RGEs of the 
modified problem in order to do so): we relax one final, say, boundary condition, find the 
resulting flow for many initial values of a coupling that is not fixed initially, and finally 
reject flows that do not satisfy the final boundary condition of the original problem. If 
multiple flows remain, then the original problem has multiple solutions. 

It is a much harder problem, in general, to count the total number of solutions for a 
given system of RGEs and boundary conditions. It can be done, in principle, by relaxing 
all of the final boundary conditions and instead scanning over all initial possible initial 
values of parameters that are not fixed by the initial boundary conditions of the original 
problem. In this way, we convert the original problem to an initial value problem, for which 
uniqueness of the flow is guaranteed (assuming the equations are sufficiently well-behaved) 
by a theorem. 

While this is certainly easy enough for a 2-parameter flow as above, it is out of the 
question for the CMSSM in which, as we will see, we would have to scan over 8 GUT-scale 
MSSM parameters for a single point in CMSSM parameter space. Thus we cannot be 
sure of finding all solutions, and therefore it is essentially impossible for us to conclusively 
declare that a generic point in the CMSSM parameter space is ruled out, in the absence of 
a theorem on the number of solutions. 

The CMSSM is, moreover, further complicated by three features. Firstly, there are 
really three 'endpoints' to the flow, namely the GUT, SUSY, and weak scales. Secondly, 
some of the boundary conditions simultaneously specify the locations of the endpoints 
themselves, as well as the values that the Lagrangian parameters take at those endpoints. 
Thirdly, the boundary conditions are themselves non-linear relations among the Lagrangian 
parameters, meaning that multiple solutions could arise even if the RGEs were linear.^ 
Nevertheless, we can identify multiple solutions easily enough in the following way. Given 
an algorithm that finds one solution, relax one of the boundary conditions (as we did for 
the BKT example) and scan over one of the couplings appearing in it, finding one solution 
for each point. If more than one of these points satisfies all of the original BCs, then the 



^In the BKT case, where the RGEs are non-hnear, but the boundary conditions (BCs) are hnear, the 
existence of multiple solutions can be attributed unambiguously to the RGEs. Similarly, if the RGEs are 
linear and the BCs are non-linear, multiple solutions can be said to arise from the BCs. But in cases like 
the CMSSM, where both RGEs and BCs are non-linear, there is no sense in which multiple solutions can 
be attributed to one or the other. 



original problem has multiple solutions. 

3 Multiple solutions in the CMSSM 

The CMSSM [9-13] is a softly broken model of supersymmetric field theory, with bound- 
ary conditions at three scales and is therefore subject to the possibility of several discrete 
solutions, even though the number of boundary conditions plus the number of input pa- 
rameters is equal to the number of free parameters. It remains the phenomenologically 
most studied example of assumptions about supersymmetry breaking terms in the MSSM, 
and is still of high interest. In this section, we will find some of its multiple solutions, map 
out the regions of parameter space in which we can find multiple solutions, and illustrate 
their properties. However, as we have already discussed, we can by no means guarantee 
that we have found all of its solutions. 

The CMSSM has the following parameters: tan/3, the ratio of the two Higgs vacuum 
expectation values; sign(/u), the sign of a parameter in the MSSM superpotential; mo, 
which is equal to the SUSY breaking family and flavour universal scalar mass terms in the 
Lagrangian; Mi/2) which is equal to the gaugino mass SUSY breaking mass parameters; 
and Aq, which sets the SUSY breaking trilinear scalar couplings (these are equal to the 
analogous Yukawa matrix multiplied by Aq). The soft SUSY breaking terms are nearly 
all set at the GUT scale, which is defined to be the scale where the gauge couplings are 
unified, typically ~ 2 — 3 x 10^^ GeV. Usually, one expects to obtain two solutions: one for 
each sign of fi, although occasionally one runs up against physical boundaries (such as an 
unstable electroweak minimum), and one or both of these solutions is unphysical. 

3.1 CMSSM boundary conditions 

Before giving the boundary conditions in some detail, we provide a rough sketch of them 
(and of a resulting RG fiow) in Fig. 2 for a certain CMSSM point. Black dots show the 
boundary conditions: at the low scale, we have boundary conditions on gauge and Yukawa 
couplings coming from experimental data. In the Figure, AfsusY is the scale at which 
superparticle masses are calculated and the Higgs potential is minimised. Thus, at MsusY; 
we have a boundary condition on /u(Msusy) coming from the Higgs potential minimisation. 
At the high renormalisation scale, we have a boundary condition on the SUSY breaking 
terms (we have illustrated this with the gaugino masses Mi). The high scale itself, Mgut, 
is usually set to be the scale Q at which gi{Q) = g2{Q), providing another boundary 
condition ((73(Mgut) is close to these, but not exactly unified in the MSSM; this small 
discrepancy can easily be explained by GUT threshold corrections [14, 15]). The boundary 
conditions are linked by the RG flow, a set of differential equations in these variables. 

The RG equations of the MSSM are non-linear, coupled, homogeneous first order 
equations. For example, at one-loop order we have equations governing the evolution of 
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Figure 2. Boundary conditions and renormalisation group flow at CMSSMIO.1.1 [16]. The abscissa 
is a logarithmic scale. 



the two Higgs soft SUSY breaking mass parameters tuhi ^^'^ ^-^2 of the form [17] 
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where t = InQ. There are other similar equations for each of the quantities on the right- 
hand side of Eq. (3.2). 

We specify our boundary conditions explicitly below for: the Yukawa couplings of the 
top, bottom, and tau, viz. ht, hh, and hr, respectively; the 3 MSSM gauge couplings, 
gi, i E {1,2,3}; the various SUSY breaking scalar masses, m^; the gaugino masses, Mj, 
i G {1, 2, 3}; and for the SUSY breaking trilinear scalar couplings for top, bottom, and tau. 
At, Ah and Ar- Other boundary conditions include those on the parameter fi appearing in 
the superpotential^ W D 11H1H2, and on ?tt-|, that mixes the two Higgs doublets in the 
potential V D m'^H2Hi. These last two boundary conditions come from the minimisation 
of the Higgs potential with respect to the neutral components of Hi and H2. 



*The circumflex indicates a superfield. 



In all, we have the boundary conditions 

tan l3{Mz) = tan/3(input) (3.3) 

ht[Mzj - , . . , hh,r[Mz) - , . ^, (3.4j 

v{Mz)s\nf3 v[Mz) cos p 

viMz) = ^,^^mB£L= (3.5) 

{^lgfiMz) + 9i{Mz) 

gi{Mz) = gi{ex.p), 52(Afz) = 5'2(exp), 5r3(Mz) = 53(exp). (3.6) 

MsvsY = ym^^iMsvsY)'mf^iMsvsY) (3.7) 

/i^(MsusY) = 2 (tan2/3(MsusY) [m2^^(MsusY) tan/3(MsusY)- 

ml^iMsvsv) cot/3(MsusY)] - M|(Msusy)) , (3.8) 

mi(MsusY) = "^"^^(^susy) (^|^(MsusY)+m|^(MsusY) + 2^^(MsusY)) (3.9) 

9i{Mgvt) = 92{Mgvt) (3.10) 

Mi(Mgut) = M2(Mgut) = M3(Mgut) = Ml/2 (3.11) 

m|(MGUT) = mJ(MGUT) = rnKMcux) = m|(MGUT) = m|(MGUT) = mg/3 (3.12) 

= m^^(MGUT) = m^^(MGUT) = ^o (3-13) 

^«(Mgut) = ^0/3, ^j(Mgut) = ^0/3, ^e(MGUT) = A0I3. (3.14) 

The parameters in Eqs. (3.3)-(3.14) are in the modified dimensional reduction (DRED) 
scheme [18]. We have labelled the input parameter tan/3 as tan/3(input). The parameters 
mb^t,T{Mz),gi,2,3iMz), and Mz{Mz) are obtained from experimental data, subtracting 
loops due to sparticles and Standard Model particles. The Standard Model electroweak 
gauge couplings gi(exp) and g2{ex.p) are fixed by the fine structure constant a and the 
Fermi constant, Gp. The values of 5fi(exp) are corrected by one-loop corrections involving 
sparticles. The parameter v{Mz) ~ 246 GeV denotes y^j^ipW^y+^ufCM^), where vi and 
V2 are the vacuum expectation values of the neutral components of the Higgs doublets Hi 
and H2, respectively. The '(exp)' denotes that the value derives from experimental data. 
m?> = m^. — ti/vi are fixed by the soft SUSY breaking mass parameters for the Higgs 
fields m'jj, i G {1,2}, as well as by the tadpole contributions ti coming from loops. These 
tadpole contributions have terms linear in /x(MgusY) as well as terms that are logarithmic 
in it, and so Eq. (3.8) is not a simple quadratic equation for /i(MgusY)- -^3 is a 3x3 matrix 
in family space. For further details, see the SOFTSUSY manual [19]. 

Spectrum calculators for the MSSM that are currently in the public domain, namely 
ISASUSY [20], SOFTSUSY [19], SPheno [21], and SUSPECT [22], solve Eqs. (3.3)-(3.14) (or 
equations very similar to them) and thus find 'the' RG flow by an iterative method. The 
particular algorithm used by SOFTSUSY, for example, is shown in Fig. 3.1. Iteration can 
only find at most one solution, for a given starting point: guesses for MSSM parameters 
are initially input to step (a), and the algorithm is applied. If a self-consistent solution to 
the whole system of boundary conditions and RG fiow is found, the parameters at step (c) 
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Figure 3. Iterative algorithm used to calculate the SUSY spectrum. Each step (represented by a 
box) is described in more detail in the text. The initial step is the uppermost one. 



remain approximately the same on successive iterations, and the iterative loop is exited, 
returning a single solution. If we are to find multiple solutions, this iterative algorithm 
must be modified. 

3.2 Finding some multiple solutions 

In order to exhibit multiple solutions, we follow the aforementioned prescription of changing 
the boundary conditions slightly. We leave all of the boundary conditions described above 
unchanged, except for the one for /z(-/VfsusY), which we allow to float instead. Through 
the Higgs minimisation condition Eq. (3.8), this allows a prediction M^(pred) for the Z^ 
boson pole mass: 



M| (Msusy) 



M| (pred) 



m 



■H2 



tan2/3(MsusY 
-2fi^{MsvsY), 
Mz(Msusy) + hI^(Msusy), 



(Msusy) tan ^(Msusy) 



- m 



5,(MsusY)cot/3(MsusY 

(3.15) 
(3.16) 



where H^^(Msusy) are the radiative corrections to the Z mass which can be found in 
Ref. [23]. When Eq. (3.16) agrees with the experimentally determined central value, 
Mz{exp) = 91.1887 GeV, we have a consistent solution of the system of boundary condi- 
tions and renormalisation group equations. Thus, in the algorithm, we supplant Eq. (3.8) 
by Eqs. (3. 15), (3. 16). This is much the same approach as the one we took in the BKT toy 
example above, in that we have relaxed a boundary condition, viz. Eq. (3.8) and scanned 
over a coupling, viz. /^(Msusy)) that appeared in it. Our algorithm can still find at most 
one solution for each value of /i(MgusY), but more than one value of /i(MgusY) might 
satisfy the original boundary conditions. 
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Figure 4. Multiple solutions in the CMSSM at mo = 2800 GeV, M1/2 = 660 GeV, Aq ^ 
and tan/3 = 40. Along the curve, all boundary conditions are applied except for the Higgs 
potential minimisation condition for /i(AfsusY), Eq- (3.8). Where the curve intersects the line 
M2(pred)/M^(exp) ~ 1 (marked as A, B and C), there are solutions that are consistent with all 
of the boundary conditions. The insert shows a zoomed region of the plot. 



We thus use a modified version of S0FTSUSY3 . 3 . 7 with the altered algorithm, and 
Standard Model input parameters as listed in Appendix A. We show M^(pred) as a function 
of /x(Msusy) for one such CMSSM point in Fig. 4. For ^ < —550 GeV, there is no physical 
solution due to ^0 being tachyonic (signalling that the desired electroweak minimum is not 
present). In the Figure, we see the usual two solutions (one for /i > at point B and one 
for ^ < at point C, both of which S0FTSUSY3 . 3 . 7 finds) plus an additional solution (at 
the point A) for /x < 0. Along the curve, Mx, gauge and Yukawa couplings, and Higgs 
soft mass parameters vary, indeed all MSSM parameters that are not boundary conditions 
vary. Between /u(MgusY) = —548 and /u(Msusy) = —530, the first term on the right- 
hand side of Eq. (3.15) is monotonic, and so it is only its combination with -2/i^(MsusY) 
that renders the minimum. We have strong evidence to suggest that there are no other 
solutions at larger values of |/i(MsusY)h there, the second term in Eq. (3.15) dominates, 
and M|(pred) is negative. 

The multiple solutions are not due to non-linearities in Eq. (3.15) alone. We have 
shown this by taking each of our solutions in turn (say. A) and then scanning over /^(Msusy) 
while not performing the RG flow, or applying the other boundary conditions, calculating 
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quantity 


solution A 


solution B 


solution C 


M^o/GeV 


282 


282 


281 


M^o/GeV 


502 


497 


471 


M^o/GeV 


558 


548 


510 


M^o/GeV 


610 


605 


593 


M +/GeV 


503 


497 


470 


M+/GeV 

A,2 


609 


604 


592 


mg/GeV 


1612 


1612 


1612 


/i(MsusY)/GeV 


-545 


-535 


497 


"ii(MsusY)/105 GeV2 


0.800 


0.809 


1.07 


m2,^(MsusY)/105 GeV2 


-1.94 


-1.83 


-1.42 


ht{MsvSY) 


0.840 


0.839 


0.836 


At{MsvsY)/GeY 


-1056 


-1057 


-1064 


Mx/lQis GeV 


1.94 


1.93 


1.89 


51 (Mz) 


0.460 


0.470 


0.456 


g2{Mz) 


0.634 


0.640 


0.633 



Table 1. Differences in CMSSM parameters and spectra for the multiple solutions of the parameter 
point mo — 2.8 TeV, M1/2 = 660 GeV, tan/3 = 40 and Aq — displayed in Fig. 4. We display some 
masses and parameters of interest for the different solutions. Above the central horizontal line, all 
masses are pole masses, whereas below the line, all quantities are evaluated in the DRED scheme. 
Out-of-the-box S0FTSUSY3 . 3 . 7 finds solutions A and C. 



M^(pred)/M^(exp). This system is observed to have only one solution for a given sign of 

Fig. 4 displays some non-smooth kinks at // ~ ±50 GeV: these occur when |M o| + 
|M o| = Mz- For |/i| > 50 GeV, a non-smooth piece is introduced in the one-loop self- 
energy of the Z^, II^^{Ms\jsy) as the particles in the loop are no longer on-shell. The 
inverted spike at /x(Msusy) ~ 6 GeV can lead to multiple solutions for other values of mo, 
-^1/2 1 ^0 and tan/3. Here, the charginos are approaching zero mass, changing the value 
of the running values of the MSSM electromagnetic coupling as extracted from data. This 
changes the electroweak couplings, which changes Mgut, significantly changing in turn 



m 



^^(Msusy) in Eq. (3.15), which dominates the prediction for M2(pred). 



We display the respective spectra of solutions A,B and C in Table 1. The spectra 
show some notable differences, illustrating the fact that the solutions are physically dif- 
ferent, leading to the possibility of their discrimination by collider measurements. Masses 
whose tree-level values depend upon the value of ^, such as the heavier neutralino and 
chargino masses, show the largest differences. Other sparticle and Higgs masses do have 
small per-mille level differences. We also see some differences in the DRED scheme running 
parameters between the solutions. We have found other points in CMSSM parameter space 
with several solutions where some sparticle masses differ by hundreds of GeV, but the ad- 
ditional solutions had sparticles lighter than Mz/2 and were obviously phenomenologically 
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Figure 5. Multiple branches of solutions in the CMSSM at M1/2 = 660 GeV, Aq = GeV and 
tan/3 = 40. Each point on the curve is a solution consistent with all boundary conditions. The 
value of Too ranges as per the horizontal axis, and corresponds to a horizontal slice through Fig. 7. 
The curve exhibits a single solution for to-q < 2796.7 GeV and multiple solutions (3 then 2 then 4) 
as rriQ increases above this value. 



excluded by LEP, which saw no significant evidence for sparticles in millions of Z^ decays. 

The CMSSM parameter point at which we have found our multiple solutions A, B, C 
is by no means alone. In Fig. 5, we scan in mo as well as /i(MsusY) in order to show the 
appearance and disappearance of various branches of the solutions. We see the appearance 
of multiple solutions for m-o > 2796.7 GeV, which then collapse to the usual two solutions, 
one for ^ > and one for fi < 0, when ttiq = 2810 - 3400 GeV. For mo = 3400 - 3700 
GeV though, there are four solutions. In Fig. 6 we exhibit another point in parameter 
space that has multiple branches of solutions which predict a CP-even lightest Higgs mass 
consistent with recent LHC measurements of a Higgs boson [24, 25]. For mo = 7572 — 7595 
GeV^ there are 3 solutions, 2 with ^(Msusy) < 0, which have m^o in the range 125.4 GeV 
to 125.7 GeV. 

We now wish to investigate how prevalent in CMSSM parameter space the multiple 
solutions are. In Fig. 7, we fix tan/3 and Ao, and scan over m,o and M1/2, showing by 
colour how many solutions (with either sign of /u) we are able to find at each point. White 
regions correspond to no solution (because, for example, there is no electroweak symme- 
try breaking), and red, orange, yellow and green correspond to 1, 2, 3, and 4 solutions, 
respectively. Grey regions have not been scanned, but the large grey area on the right 

''As is well known, getting nih ~ 125 GeV in the MSSM in general requires unnaturally large stop quark 
masses, and hence large mo in the CMSSM. 
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Mi/2=1 TeV, tan|3=40, Ao=-5 TeV 
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Figure 6. Multiple branches of solutions in the CMSSM at M1/2 = 1 TeV, Aq = -5 TeV and 
tan/3 — 40. The value of mo ranges as per the horizontal axis. We plot here the predicted values 
of the Higgs mass, m^a . There are 3 solutions in the range 7572 GeV < itiq < 7595 GeV which are 
consistent with recent LHC measurements of a Higgs boson mass. 



of the plot is in the usual region of unsuccessful electroweak symmetry breaking, and we 
expect solutions. The large red region has a single solution for /i > 0; in this region 
there is no solution for ^ < due to tachyonic Aq. The orange region has the usual 2 
solutions, one for each sign of fj,. Beyond this, near the boundary of electroweak symmetry 
breaking, multiple solutions for both signs of /x are prevalent. The solutions found near 
the edge of successful electroweak symmetry breaking correspond to the spiky apex of the 
curve in Fig. 4 shifting down to intersect the horizontal dotted line in 4 distinct places, 
corresponding to 4 correct predictions for Mz- 

Fig. 5 shows the value of // for a slice through Fig. 7 at M]^ ^ = 660 GeV. For low 
mo we are in the red region and there is only 1 solution. As we increase mo, just before 
entering the region with 2 solutions, we hit a small patch of parameter space which exhibits 
3 solutions, 2 of which are for /x < 0. The points in parameter space where this occurs 
are indicated (for those with keen eyesight) in Fig. 7 by the isolated yellow dots between 
the red and orange regions. We have verified that these regions of 3 multiple solutions are 
not confined to single points in parameter space, but rather occupy a finite area, as shown 
by the inset in Fig. 5, and as shown by the small region ttiq = 7572 — 7595 in Fig. 6. It 
is expected that this area runs the length of the boundary between the red and orange 
regions. Whilst most of the multiple solutions are fairly near the boundary of electroweak 
symmetry breaking, there is nevertheless a significant volume of parameter space where 
the multiple solutions play a role, and they should not be ignored in phenomenological 
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Figure 7. Parameter scan of the CMSSM, with colour indicating number of solutions at a given 
point: white, red, orange, yellow and green being 0, 1, 2, 3 and 4 solutions respectively. The red 
region at the top-left has 1 solution for /x(Msusy) > and no solutions for ^(A^susy) < 0. The 
orange region to the right of this red region has 2 solutions, 1 for each sign of /J.(AfsusY)- Further to 
the right are regions with 3 and 4 solutions. Phenomenological constraints have not been applied. 



analyses. 

We have performed further scans of parameter space at selected values of tan (3 and Aq, 
and counted solutions in the mo-Mi/2 plane. The qualitative features of these other slices 
through the parameter space of the CMSSM are similar to Fig. 7, including the existence 
of a region with three solutions between the red and orange regions. One occasionally 
finds small areas, close to the edge of successful electroweak symmetry breaking, that have 
six solutions, two for /i < and four for ^ > 0. For example, tan/3 = 20, ^o = —1000 
GeV, Ml/2 = 170 GeV and ttiq = 2840 GeV is a point with six such solutions. The 
regions with three and four solutions occur at larger values of mo close to the region 
of electroweak symmetry breaking. While the position of the no electroweak symmetry 
breaking boundary is extremely sensitive to the top quark mass mt, so that uncertainties 
on it remain high [26, 27], we have checked that the additional solutions remain (but move 
in mo) when we change m^t- 

4 Discussion 

We have shown that quantum field theories with multiple parameters that have bound- 
ary conditions imposed at multiple renormalisation scales may admit several physically 
distinct solutions for the same boundary conditions. Each separate solution corresponds 
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to a different RG trajectory, consistent with the RG flow itseff being unique. We have 
exemphfied this analyticahy with the simple BKT model, and shown the relevance of mul- 
tiple solutions to high energy physics, where each solution generally corresponds to distinct 
phenomenology. 

Imposing boundary conditions on the CMSSM naturally involves different renormal- 
isation scales, namely the weak scale, MgusY; and the gauge coupling unification scale. 
Thus, the CMSSM is a good candidate for a theory possessing several different solutions, 
for a single set of input parameters sign(;^), ttiq, Mi/2) ^o and tan/3. Publicly available 
computer programs solve for the RG flow using iteration, which flnds at most one solution, 
given the boundary conditions, for a given point in the CMSSM parameter space. How- 
ever, we have shown that multiple, physically distinct solutions exist in extended regions 
of parameter space. We showed this by relaxing a boundary condition, scanning over the 
parameter ^(Msusy), and (for reasons of computational efficiency) completing the rest 
of the calculation by iteration. We find additional, previously unknown and unexplored 
solutions at different values of fJ-iMsusv), which we cryptically referred to in the title as 
"the dark side of the //" . However, we might well have found yet more solutions if we had 
scanned over other parameters, so in reality all of the parameters that are not directly fixed 
by a BC have just as much of a dark side, waiting to be explored. 

In order to be sure of finding all of the solutions, one could turn the boundary value 
problem into an initial value problem and perform the RG flow in one direction only. One 
possible method would be to start at the GUT scale and, for a given sign(/i), niQ, Aq, M^^, 
and tan/3(input), use Mgut, 5i(Mgut) = 52(Mgut), 53(Mgut), /li(MGUT), /ib(MGUT), 
/it-(Mgut)) A*(-Mgut) and m^MGUT) as scanning parameters. This 8-dimensional scan, if 
performed flnely enough, would then find all of the possible solutions that match 5'i(exp), 
g'2(exp), g'3(exp), ruh, rrit, m-r, Mz and tan f3{Mz)-^ Note that However, we found that even 
while scanning only in one dimension and finding the rest of the parameters by iteration, 
some of the low-energy output parameters are rapidly varying functions of the high-energy 
inputs, so that finding some of the solutions is computationally intensive. This problem 
would inevitably be more acute for an 8-dimensional scan, and we conclude that it would 
be extremely difficult to know for sure whether one has found all of the solutions, in the 
absence of a theorem as to their number. 

It seems at least plausible, then, that some regions of CMSSM parameter space have 
been ruled out erroneously, because of the existence of additional solutions, yet to be found, 
with a phenomenology that is consistent with current bounds. For example, we have seen 
that neutralinos and chargino masses may differ at the tree-level (i.e. at order 100% level), 
between the multiple solutions, and so bounds coming from dark matter searches are liable 
to change, since they are very sensitive to the mass of the dark matter candidate (in this 
case, the lightest neutralino). Similarly, we have seen in Fig. 6 that the same parame- 
ter point can lead to multiple predictions for m^o; in the coming era of precision Higgs 
physics, we may well find points for which the solution found by the standard algorithms 



®tan/3(Mz) is, strictly speaking, an output of our scan. Nevertheless, we follow convention and include 
it in the above list of input parameters for the CMSSM. 
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is incompatible with the measured m/,, while other solutions are compatible. 

Should we therefore consider all existing exclusions of the CMSSM parameter space as 
suspect? Fortunately, there are sectors in which we do not expect the phenomenology to 
be greatly changed for different solutions, such that existing exclusions from experimental 
searches may be considered robust. Important examples are collider searches for gluinos 
and squarks decaying into jets and missing energy, which provide powerful constraints 
upon the CMSSM [28, 29]. The gluino/squark production cross-sections relevant for the 
calculation of these bounds depend upon the gluino and first two-family squark masses. 
Since these are fixed at the GUT scale by the boundary condition M3(Mgut) = -^1/2 and 
mq;(AfGUT) = iTT'O, they only differ between the different solutions at the loop level. Pertur- 
bativity implies that these corrections are small, and so one expects that the production 
cross-sections will not differ greatly between the various solutions. Thus, we expect the 
previously calculated limits on the CMSSM from squark and gluino production to approx- 
imately hold. One significant caveat to this is in regions of parameter space where the 
ratio of the neutralino mass to the gluino or squark mass differs significantly between solu- 
tions. This can significantly change the kinematics of the decay of the gluino and squark, 
and therefore can affect the acceptance of signal events passing the cuts applied in the 
experimental analyses. 

For the time being, in the absence of a signal for supersymmetry at colliders, multiple 
solutions add a potential loop hole to some of the exclusion bounds derived in models such 
as the CMSSM. If evidence for supersymmetry is found in the future, the phenomenon will 
present a new facet: points in the parameter space of candidate SUSY theories might be 
incorrectly discarded on the grounds of not being able to explain the signal, when other, 
unknown solutions, corresponding to the same parameter point, can. We thus might be led 
in the wrong direction in searching for the correct theoretical description of new physics. 
Conversely, if all the solutions are known, then we should be able to discriminate between 
them given sufficient measurements, given that they are physically distinct. 

There has been a recent growing industry in multi-dimensional phenomenological fits 
of parameter space in the CMSSM to collider and astrophysical data (see Refs. [30-35] 
for some examples). These fits use either Bayesian or frequentist statistics. The Bayesian 
fits yield posterior probability densities, which are weighted by the probability masses in 
marginalised parameters. The addition of multiple CMSSM solutions to these fits would 
only have a negligible affect if the probability mass associated with the non-standard so- 
lutions is negligible. Performing a Bayesian fit with the additional 8 dimensions of our 
scan would be an interesting exercise, to see what the effect of the multiple solutions is. 
Such fits, since they scan in m^ and /i (among other parameters) rather than tan/3 (and 
fixing ^ by Mz {ey^p)), will implicitly and automatically take into account a Jacobian factor 
that was calculated for this purpose [36, 37]. Frequentist fits must scan in the additional 
solution space and may be even more vulnerable to changes if they have a better best-fit 
point, since parameter estimation always relies on A^^, the difference in x^ between the 
current parameter space point and the best-fit point. It remains to be seen how much these 
effects will change the fits in each case, but at the moment they remain incomplete without 
the inclusion of the multiple solutions exhibited here. 
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Finally, we emphasise that the CMSSM is only one example of a model where the 
SUSY breaking conditions are mostly set at a scale much higher than the weak scale; we 
expect these other theories to feature multiple solutions as well. 
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A Standard Model input parameters 

As input parameters, we have [38] Gp = 1.16637 GeV~^, pole masses mt = 173.5 GeV 
and rrir = 1.77669 GeV, modified minimal subtraction scheme Standard Model values of: 
rribimb) = 4.18 GeV, Mz(exp) = 91.1887 GeV, a{Mz)~^ = 127.916, asiMz) = 0.1187. 
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